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ABSTRACT 

The greatest challenge in the interpretation of galaxy clustering data from any surveys is galaxy 
bias. Using a simple Fisher matrix analysis, we show that the bispectrum provides an excellent 
determination of linear and non-linear bias parameters of intermediate and high-z galaxies, when all 
measurable triangle configurations down to mildly non-linear scales, where perturbation theory is still 
valid, are included. The bispectrum is also a powerful probe of primordial non-Gaussianity. The 
planned galaxy surveys at z > 2 should yield constraints on non-Gaussian parameters, and /^^, 
that are comparable to, or even better than, those from CMB experiments. We study how these 
constraints improve with volume, redshift range, as well as the number density of galaxies. Finally, 
we show that a halo occupation distribution may be used to improve these constraints further by 
lifting degeneracies between gravity, bias, and primordial non-Gaussianity. 

Subject headings: cosmology: theory - large-scale structure of the Universe 



; 1. INTRODUCTION 

Why study high-z galaxy surveys? The recognition 
that bary on acous tic oscill ations in the galaxy power 
, spectrum (ICole et a l. 2005: Eisenstein et alJ[2005t iHtitsl 
■ 120061 : iPercival et al.ll2007l ) are an excellent probe of the 
\ nature of dark energy has led to several proposals for 
large- volume redshift surveys at z > 1. 

Galaxies are a biased tracer of the underlying matter 
distribution. The use of highly biased tracers, such as lu- 
minous red galaxies at lower z and Lyman break galaxies 
or Lyman-a emitters at higher z requires a reliable mod- 
\ elling of non-linearity and scale-dependen ce of galaxy 
. bias, even at relatively large spatial scales (jSmith et al.l 
; [20M [20071 : rMcDonaldl [2006). 

The distribution of galaxies is non- Gaussian. The 
galaxy bispectrum, the three-point correlation function 
in Fourier space, does not vanish. It has been known 
for more than a decade that the bispectrum is an excel- 
lent tool for measuring galaxy bias parameters, indepen- 
dent of the overall normalization of dark matter fluctua- 
\ tionsJFry 1994; Matarrcsc et al. 1997; Scoccimarro et ahl 
< l2001al ). This method has been applied successfully 
to existing galaxy surveys such as the 2dFGRS and 
the SDSS, y ielding cons t raints on non-linearities in 
galaxy bias fVerde et al.' '2002- 'Pan fc SzapTidil l2005l : 
iGa ztahaga et al. 2005; Nishimichi et al. 2006|) as well 
as on the Halo Occupation Distribution (HOD) 
(jKulkarni et al.l[2007l ). Moreover, it has been shown that 
the galaxy bispectrum contains additional cosmological 
information that is not present in the p ower spectrum 
(jSefusatti &: Scoccimarr o 2005; Sefusatti et al. 2006). 

The galaxy bispectrum on large scales, or other sta- 
tistical tools that are sensitive to the higher-order cor- 
relation of galaxies, are sensitive to statistical properties 
of primordial fluctuations: primordial non-Gaussianity 
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Scoccimarro ct al. 2004; Hikagc ct al. 2003). The com- 
mon belief is that the cosmic microwave background 
(CMB) is most sensitive to primordial non-Gaussianity; 
however, the number of Fourier modes available in the 
galaxy survey data is much larger than that in the CMB 
data, as the former probes the full three-dimensional 
structure of density fields. The galaxy bispectrum can, 
therefore, outperform the CMB bispectrum, provided 
that the other sources of non-Gaussianity, such as non- 
linear gravitational evolution and galaxy bias, are under 
control. Furthermore, it should always be emphasized 
that galaxy surveys provide information on the spatial 
scales that are much smaller than those probed by CMB; 
thus, these two measurements are complementary to each 
other. 

Motivated by these considerations, in this paper we 
study how well the planned high-z galaxy surveys would 
constrain primordial non-Gaussianity and galaxy bias 
parameters. High-z surveys are more useful for this task 
than low-z ones owing to much weaker non-linearities 
in matter clustering and redshift-space distortion, which 
allows us to use the galaxy bispectrum to the smaller 
spatial scales that are inaccessible by low-z surveys. 

This paper is organized as follows. In section [2] we 
describe the leading contributions to the bispectrum of 
the observed galaxy distribution. In section [3] we present 
our Fisher matrix analysis for the galaxy bispectrum, and 
in section [4] we show our predictions for constraints on 
galaxy bias and primordial non-Gaussianity from sample 
galaxy redshift survey designs, studying dependence on 
the survey volume, maximum wavenumber, number den- 
sity of galaxies, and redshifts. We also present our pre- 
dictions for the ongoing, upcoming and planned galaxy 
surveys. In section[5]we extend our fisher matrix analysis 
to the halo occupation distribution as a tool to describe 
galaxy biasing at large scales. Finally, in section [H] we 
present our conclusions. 

2. THE GALAXY BISPECTRUM 
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2.1. Primordial non-Gaussianity 

We explore two parametrizations of primordial non- 
Gaussianity which are motivated by inflationary models. 
While representing a wide variety of non-Gaussian mod- 
els, these parametrizations are by no means exhaustive. 
Our method can be appli ed to any other functional forms 
of th e bispectrum (e.g.. iLiguori et ahl 120061 : IChen et all 
[200l in a straightforward way. 

2.1.1. Local model 

The first one is described by the local expression for 
Bardeen's curvature per turbations during the matter era, 
$(a;) , in position space (iGangui et"aLlll994l ; IVerde et al.1 
[2000: IKomatsu fc Sperge]|l200in 



<I>(x) = cI>G(x) + /jv°i[1>^(x) 



where $g(x) is a Gaussian field and is a con- 

stant characterizing the amplitude of primordial non- 
Gaussianity. In this case, the leading contribution in 
the expansion to the bispectrum, B^^'^'^^ki , k2 , k^) , 
of the curvature field is given by 



(1) 



t2f]^l[P^{kl)P^{k2) 
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where we approximate P^{k) ~ and 



C<5 — 



(2) 



(3) 



which quantifies departure from a scale-invariant spec- 
trum. We include Hs explicitly, as we are interested 
in determining how a departure fr om scale invariance , 
which has been hinted by WMAP (jSpergel et al.ll2006f ). 
would affect detectability of primordial non-Gaussianity 
in the distribution of galaxies. For this local model, most 
of the signal is given by squeezed triangular configura- 
tions, ki <C ^2, fca. The local type of non-Gaussianity 
described by equation is predicted in models such 
as the curvaton scenario (iLvth et al.l l2G03*). models with 
inhomogeneous reheating ( Dyali et a l. 2004b *fl), multi- 
ple field inflationary models ( Bernardeau fc Uzanll200^ 
or generically in models where the non-linearities arise 
from the evolution of perturbations outside the horizon. 

The best limits to date on possible values for 
the parameter come from measurements of the 

microw ave bac kground bispectr um on the WMAP 
(IKomats u et al. 2003; iSpergel et all l2006l : 



data 



iCremineUi et~all [200l -36 < f]^l < 100 at 95% 
C.L., which corresponds to the l-cr error of 



A/]^2-=34 (WMAP3), 



(4) 



which is a factor of 50 better than the limit from COBE 
(|Komatsu et al.ll2002[) . Upon completion, WMAP is ex- 
pected to reach A/]^£- ~ 20, while the Planck satel- 



lite would yield A,/ ^ r, — y n-Lvuinciua u <x, uusu 
iBabich fc Zaldarriagall20'ol lYadav et al. 200"^. 

Measurements of the galaxy bispectrum in the SDSS 
main sample are expec ted to yield Af^Sf. — 150 
Scoccimarro et aTl l2004[ ). Recently iPillepich et al.l 



floe. 



3 (IKomats u fc Sp ergell I2OOII : 



2006[ ) pointed out that a full-sky measurement of the 



bispectrum of fluctuations in the 21-cm background 



might reach / l yr,' — 1 , while a more aggressive analy- 
sis by ICooravl ()2006f ) shows that the same observations 
could reach A/]^£' ~ 0.01 in principle. Other large- 
scale structure probes such as cluster abundance, on the 
other hand, is unlike ly to improve CMB limits on 
(jSefusatti et al.l[2007l ): however, it should provide an im- 
portant cross-check of the results if a signiflcant 
was detected in the CMB, and it should not be forgotten 
that the spatial scales probed by the cluster abundance 
is smaller than those probed by the CMB. 

2.1.2. Equilateral model 

The second model for primordial non-Gaussianity is 
given by 
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IBabich et all (|2004D and iCreminefli et all (|2007D have 
shown that this form provides a good approximation to 
the bispectra predict ed by higher der i vatives and DBI 
inflat ionary models (jCreminellii l2003t lAlishahiha et aTl 
I2OOI . 

The bispectrum in equation ([5]) is normalized in a such 
a way that for equilateral configurations {ki = k2 — k^ — 
fc), it coincides with the local form given in equation ([2]). 
The important difference is that this form has the largest 
contribution from the equilateral configurations, as op- 
posed to the local form in which the largest contribution 
comes from the squeezed configurations. The current 
limits from WMAP ar e -256 < f^^ < 332 at 95% C.L. 
(jCreminelli et al.ll2007| ). which corresponds to the l-cr er- 
ror of 

A/^\ = 147 (WMAP3). (6) 



2.1.3. The primordial density bispectrum 

Density fluctuations in Fourier space, ^k, are related to 
the curvature perturbations, <i>k, via the Poisson equa- 
tion, (5k(a) = M(fc;a)$k, where 



M{k,a) = - 



_ 2_D(a) , 2 
3Hm 



k^T{k). 



(7) 



Here a is the scale factor, T(fc) is the matter transfer 
function, and D(a) is the growth function ^ . 

This allows us to write the primordial contribution to 
a generic n-point function of the matter density fields in 
terms of the respective correlator of the curvature per- 
turbations as 

{6]^^5]^^...5]^^)i^M{ki\a)M(k2\a)...M{kN\a) 

x($ki$k....$k„). (8) 

^ The function , M{k;a), uses the same definition as in 
IVerde et all 112003) for Vim = 1, but it differs from M{k) given 



in IHikageetaLT i 20061 ) , where the dependence on the growth 
function has been taken out as M( k;a) = M^'^^ (k)Dia) . 
Also, the same function is defined in IScoccimarro et al.l (12004 ) 
in terms of the gravitational potential (i.e., the 0-0 component 
of the metric perturbations) during radiation domination, so that 
M{k;a) = — -^M^^^ [k; a) leading to a different definition of the 
non-Gaussian parameter, /jvl = ~^fp/if 
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In particular, the initial (primordial) matter bispectrum, 
Bi{ki,k2, fcs), is given by 

Bi{h,k2, ks) = M(fci)M(fc2)M(fc3) 

xB^{ki,k2,k3), (9) 

where we have omitted for brevity the explicit depen- 
dence of M{k; a) on a. We can also relate the linear 
density power spectrum, PL{k), to the curvature power 
spectrum, P.j,{k), as 

PL{k)^M^{k)P^{k). (10) 

A hierarchical relation between the scale dependence 
of the initial bispectrum and the power spectrum such as 
B^{k) ~ fNLP^{k) with a constant /jvl is by no means 
generic or universal. A /tvl with a peculiar scale depen- 
dence unrelated to the power spectrum appears, for in- 
stance, in string-motivated models s uch as DBI inflation 
(jAHshahiha et all l2004t IChenI l2005l ). Our analysis can 
be applied to any models of primordial non-Gaussianity, 
provided that the bispectrum can be calculated from 
those models. 

Note that a post-Newtonian effect can yield an ad- 
ditional contribution to non-linearity of primordi al per- 
turba tions and hence to non-Gaussianity (Bartol o et al.l 
I2005D . Although we do not include this effect our anal- 
ysis, it would be interesting to study how important the 
post-Newtonian effect would be for the future galaxy sur- 
veys. 

2.2. Non-Gaussianity from non-linear gravitational 
evolution 

Even if the initial perturbations are Gaussian, the sub- 
sequent gravitational evolution makes the evolved den- 
sity fields non-Gaussian. On large scales one can study 
the non-linear evolution of matter density fluctuations 
by means of perturbation theory, and write the solution 
up to the second order in 6 as 

4 ^ 4'^-f^d3gid3g25c(k-qi2)^^2(qi,q2)<5^V'5qV' (H) 

where 5*^^^ is the linear solution, and F2(ki,k2) is a 
known mathematical function given by 

^ „ , s 5 X f k] fco \ 2 n , , 

F2i).,M)^-, + -^[-i-^+f;)^-,^\ (12) 

with x = ki • k2. Therefore, one obtains 

SG(fcl,fc2,fc3) =2F2(ki,k2)PL(fci)Pi(fc2)-HCyC. (13) 

The bispectrum of matter density fluctuations (i.e., no 
galaxies yet) evolved from non- Gaussian primordial fluc- 
tuations on large scales is thus given by the sum of equa- 
tion (P)) and 

S(fci, fc2, fcs) = S/(fci, fc2, A:3) + Baiki, fc2, A^g). (14) 

As usual, we shall focus on the reduced bispectrum, 
defined as 

r,(i. I. I. \- B{ki,k2,k3) 

g fci,fc2,fc3 = ^TTrT^TTrT— — (15 

P{ki)P[k2) + eye. 

which has an advantage of being only mildly sensitive 
to cosmological parameters. That is to say, the depen- 
dence on cosmology has been "factored out" by a prod- 
uct of the power spectra in the denominator and the Qq 



component is particularly insensitive to the amplitude of 
matter fluctuations (e.g., cts)- The reduced bispectrum 
of matter density fluctuations is also given by the sum of 
two contributions: 

Q{ki,k2, fcs) =Q/(fci, ^2, ^3) + Qciki, k2, k^) 
_ Bi{ki,k2,kz) 
~ P{ki)P{k2)+cyc. 
BG{ki,k2,k3) 
"^P(fci)P(fc2)+cyc. ^ 

It is important to remember that, in the leading order, 
Qg does not depend on the linear growth factor, D{a), 
and thus it is independent of redshifts. In other words, 
Bg is proportional to D"^, which cancels D"^ in [P{k)]'^ 
in the denominator exactly. On the other hand, Qi is 
proportional to 1/D{a) because Bj (x D^, and thus it 
is larger at higher redshifts. Therefore, high-z galaxy 
surveys are expected to be more sensitive to primor- 
dial non-Gaussianity, relative to the gravitational bis- 
pectrum, than low-z ones. 

In Figure[l]we plot the equilateral configurations of the 
reduced bispectrum, Q{k) = Q(k,k,k), from non-linear 
gravitational evolution and non-Gaussian initial condi- 
tions at z = 0, 1, and 4. As mentioned earlier the local 
and equilateral model of primordial non-Gaussianity give 
the same results for these configurations. For Gaussian 
initial fluctuations, Jnl = 0, Q{k) ~ 0.57 at tree-level 
in perturbation theory and is independent of scales."* 
On the other hand, Q{k) exhibits a clear scale depen- 
dence for /atl ^ 0. a positive Jnl enhances Q{k) at 
large scales, i.e., Q{k) > 0.57, whereas a negative Jnl 
suppresses it. This is because a positive /nl results 
in positively skewed density fluctuations. It is easy to 
show that the scale dependence of the primordial com- 
ponent Qi{k) is given by 1/M(fc), where M{k) is given 
by equation As 1/M{k) oc fc~^ on large scales, we 
find that the primordial non-Gaussian signal is larger on 
large scales. This property makes it easier to find primor- 
dial non-Gaussianity in CMB observations; however, as 
there are much more modes available on smaller scales, 
the cumulative signal to noise for higher-order correlation 
functions increases due to the large number of observable 
configurations on small scales (Sefusatti & Scoccimarr^ 
2005). 

How about other configurations? In Figure [2] we 
plot the dark matter bispectrum with Gaussian or non- 
Gaussian initial conditions at different redshifts as a 
function of the angle, 9, between ki and k2, for ki — 
0.01 h Mpc~^ (top panels) and 0.02 h Mpc"^ (bottom 
panels) and k2 = 2ki. We find a marked difference in 
the configuration dependence of the two primordial bis- 
pectra, local (left panels) and equilateral (right panels), 
under consideration. 

2.3. Non-Gaussianity from galaxy bias 

The galaxy bispectrum is most useful for measuring 
galaxy bias. Assuming that galaxy formation is a local 
process and depends only on the local matter density 

* The equilateral reduced bispectrum is independent of scales 
only in the second order perturbations. A scale dependence 
arises when the higher-order terms are included (jScoccimarro et al.l 

[HID. 
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0.01 0.1 0.01 0.1 

k [h/Mpe] k [h/Mpc] 

Fig. 1. — Equilateral configurations of the reduced bispectrum of dark matter distribution at the second order ("tree-level"). The 
horizontal lines at Q{k) = 0.57 show the gravitational contribution only, which corresponds to = = The solid, long-dashed 

and short-dashed lines show the gravitational contribution with non-Gaussian initial perturbations at ^ = 0, 1, and 4, respectively. Note 
that = /jvL equilateral configurations. {Left panel) The curves above Q{k) = 0.57 show /jvl = +100, while the curves below it 

show /jv£ = —100. {Right panel) The same as the left panel but for the current WMAP3 limits on (top) and (bottom). 



field, one may expand the galaxy number overdensity, 
dg, in Taylor series of the und erlying matter overdensity, 
6, as (|Frv fc Gaztanagalll993f) 



(17) 



where bi and 62 are the linear and non-linear bias pa- 
rameter, respectively. It has been shown that this model 
describes the bispectrum or three-point correlation func- 
tions from the SDSS and 2dFGRS (e.g..lGaztafiag a et al.l 
[2b05; Nishimichi ct al. 2006), as we ll as from numerical 
simulations fe.g.. lMarfn et al.il2007[ ). The galaxy bispec- 
trum is given by 

Bg(fci,/c2,fc3)=i blB{ki,k2,k3) 

+b\b2[PL{ki)PL{k2) + eye], (18) 

up to the second order in matter density fluctuations. 
Here, B{ki, k2, fca) is the intrinsic bispectrum of the un- 
derlying matter distribution. The reduced galaxy bispec- 
trum is 

1 &2 
Qg{ki,k2,k3) ~ — (3(fci,fc2,/c3) + 72- (19) 
Oi Oi 

With Q(fci,/c2,fc3) = QG{ki,k2,k3) + fNLQi{ki,k2,k3), 
one obtains 

Qgiki, fc2, /ca) ~ — QG(fci, fe, fe) 

+ ^g/(fci,A;2,fc3) + ^, (20) 
where we have factorized f^L out from Q/ introducing 

Qi = QiUnl = 1). 

Finally, the galaxy power spectrum is given by Pg{k) ~ 
b\P{k), up to the second order in density fluctuations; 



however, corrections due to non -linear b ias appear at 
the third-order level dHe avcns ct al. logs'; iTaruval [2OOOI : 
ISmith et al.ll2006tlMcDonald.2006. ). Therefore, the bias 
parameters, 61 and &2, affect both the galaxy power spec- 
trum and bispectrum. The bispectrum helps us extract 
the cosmological information from the galaxy power spec- 
trum by providing 61 and &2- 

2.4. Redshift space distortion 

The bispectrum measured from redshift surveys is dis- 
torted along the line of sight direction by radial motion 
of galaxies. For our analysis in this paper we shall deal 
only with a spherically averaged power spectrum and bis- 
pectrum. The power spectrum in redshift space after av- 
eraging over angles in k space, Ps{k), is related to the 
real space power spectrum by 

P.(fc)=ao^(/3)Pg(fc), 

while the bispectrum is given by 

Bs{ki,k2,k3) = aQ{P)Bg{ki,k2, 

where (|Kaisedll987l: ISefusatti et"aI]|2006D 



, ^3), 



(21) 



(22) 



«o (/?) = !+ 3/3- 



^(/3) = l + -/3+-/3' 



(23) 
(24) 
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with (3 = n,n /hi. The reduced bispectrum in redshift 
space is thus given by 



Qs{ki,k2, fcs 



Km? 



Q{ki,k2,k3) + 



b\ 



(25) 



We remark that our treatment does not take into 
account a peculiar scale-dependence of redshift distor- 
tions that reduces the amplitude of non-linear corrections 
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Fig. 2. — Configuration dependence of the reduced bispectrum of dark matter distribution from Gaussian and non-Gaussian initial 
conditions, as a function of an angle, 6, between two wave vectors, ki and k2, where the magnitude satisfies k2 = 2k\ = 0.02 /iMpc~^ 
(top panels) and 0.04 hMpc~^ (bottom panels). (Left panels) = ±100. [Right panels) f'^j^ = ±200. The dotted black line shows a 

Gaussian case {/nl = 0, redshift independent), while the solid, long-dashcd and short-dashed lines show non-Gaussian cases at 2 = 0, 1 
and 4, respectively. 



to th e reduced matter bispectrum (jScoccimarro et al] 
Il999f ). This fact only partiaUy justifies the approxima- 
tion of neglecting such corrections altogether. 

3. FISHER MATRIX ANALYSIS 
3.1. Method 

In our analysis we shall consider a set of surveys char- 
acterized by their volume, V, mean galaxy density, Ug, 
and redshift range. We shall assume that these surveys 
have a simple survey geometry, i.e., a contiguous hexa- 
hedron. 

Our bispectrum estimator is given by 



(jScoccimarro et al.|[l998l ) 

B = ^ f f d\-2 f d\3 SD{qi23)S^,S^J^,, (26) 
Vb Jki Jk2 Jk3 

where the integration is over the bin defined by qi £ 
{h-Sk/2,ki + Sk/2), Vf = (27r)VV"is the volume of the 
fundamental cell in Fourier space, and 

Vb= d^qi d^q2 d^q3dD{q.i23) 

fci k2 ^3 

-Stt^ /fcifefcs Ak^, (27) 

with Afc a multiple of the fundamental frequency, kf = 
2tt/L. We assume that two coincide, i.e., AA: = kf, 
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thereby taking into account all "fundamental" triangular 
configurations. 

The variance for our estimator, to the leading order, is 
given by a triple product of the power spectra, 

Ai?2 ~ fc3 iHi P(fci)F(fc2)F(fc3), (28) 
Vb 

where S123 — 6,2,1 for equilateral, isosceles and general 
triangles, respectively. As for variance of the reduced bis- 
pectrum, we assume that variance from the bispectrum 
in the numerator dominates over that from the power 
spectra in the denominator: 

AQ^ AB^ ^^^^ 



We calculate variance of the redshift-space galaxy re- 
duced bispectrum from equations and ([^5]) with 
P{k) given by 

Pi„t(fc)^P,(fc) + -i^i, (30) 

where the second term accounts for the shot noise. We 
finally obtain 

siisk"} Ptotiki)Ptotik2)Ptotik3) 



AQiiki,k2,k3 



Vb [P.(fci)F,(fc2)+cyc.]2 



.(31) 

Once the variance of the reduced bispectrum is given, 
the Fisher matrix for a given redshift bin can be ex- 
pressed as 



dQsji) dQsji) 1 

dpa dp/3 AQ2(j)^ 



(32) 



where the parameters, Pa, represent 61, 62 and Jnl- We 
use equation (|25p for Qg, which is valid only up to the 
second order in perturbations. Neglecting higher order 
corrections would introduce systematic errors when deal- 
ing with the real data, particularly at low z and at small 
spatial scales, where non-linearity is substantial and per- 
turbation theory essentially breaks down. Since we con- 
sider high-z surveys on large scales, we expect that higher 
order effects would not affect our results very much. 

As a fiducial cosmological model we use a flat ACDM 
cosmology with matter density f2„i = 0.3, baryon density 
fib = 0.04, Hubble parameter h — 0.7, spectral index 
Us = 1 and cTg = 0.9. Of these parameters, erg and ris 
affect our forecast for the projected errors on the bias 
parameters most. We thus consider also different values 
such as those sugges ted by the WMAP 3 -yr results, ag = 
0.75 and = 0.95 (ISpergel et al.ll2006l ). 

3.2. Comments on covariance matrix 

We shall not include covariance between the cosmolog- 
ical parameters and bi, 62, and /jvi- ( We do include co - 
variance between bi, 62, and /nl-) Sefus atti et al.l ()2006f ) 
have shown that, for the SDSS main sample, an analy- 
sis with the full covariance among all parameters with a 
prior from the WMAP 3-yr results yields an error on bi 
that is twice as large as that from a simpler analysis with- 
out covariance. On the other hand, an error on 62 is not 
affected significantly. Note that the effect on 61 that they 
observed was due mainly to degeneracy between bi and 
the amplitude of matter fluctuations, as I Sefusatti et al.l 



(|2006D did not use the reduced bispectrum. We expect 
that degeneracy would be lifted in our analysis, as we use 
the reduced bispectrum in which the overall amplitude 
of matter fluctuations cancels. 

More importantly, we shall not include covariance be- 
tween different triangular conflgurations in the bispec- 
trum. The covariance arises from both observational se- 
lection functions (i.e., survey geometry and mask) and 
a connected six-point function generated by non-linear 
gravitatio nal evolution. This is a r ather crude approx- 
imation. IScoccimarro et aTl (|2004f ) have included the 
full reduced bispectrum covariance plus the peculiar sur- 
vey geometry, when they calculate the constraints on 
galaxy bias and primordial non-Gaussianity from the 
SDSS main sample. Using the same realizations of the 
survey and the same estirnator for the covariance matrix, 
ISefusatti &: Scoccimarrol ()2005D have compared the anal- 
ysis with the full covariance matrix (including the ob- 
servational selection function) and that with an approxi- 
mate diagonal Gaussian variance. They have found that 
the latter simplifled treatment overestimates the signal- 
to-noise by a factor of 2 for fcmax ~ 0.1 hMpc~^, and a 
factor of 8 for fcmax ~ 0.3/iMpc~^ at redshift zero. It 
is not clear, however, how to separate the contribution 
from non-linear evolution from the effect of the selection 
function. One would generically expect the radial contri- 
bution to be smaller at high z, as the six-point function 
from gravitational clustering becomes smaller than the 
non-connected part of six-point function (which consists 
of power spectra) at higher z. In any case, our results 
should be taken as a guide, and one needs to perform the 
full analysis including the selection functions peculiar to 
a given survey design. 

3.3. Non-linearity and maximum wavenumber 

While one can measure the galaxy power spectrum or 
bispectrum down to very small spatial scales, say, 10 kpc, 
it is challenging to extract useful cosmological informa- 
tion from such small spatial scales owing to strong non- 
linearity. Therefore, one has to decide on the maximum 
wavenumber, fcmax, below which theory may be trusted. 
Not surprisingly, since there are many more modes avail- 
able on smaller spatial scales, the amount of cosmological 
information one can extract from data grows as fcmax in- 
creases. It is therefore important to use a realistic fcmax 
in order not to overestimate the statistical power of a 
given galaxy survey design. 

How do we decide on fcmax? The first obvious thing 
to do would be to test our theory of the power spec- 
trum and bispectrum against numerical simulations. A 
value of fcmax can be found by comparing perturba- 
tion theory predictions with n umerical simulations (see 
e.g., iJeong fc Komatsul l2006l for an analysis for the 
matter power spectrum). It is likely that a simple 
model provided by the second-order (tree-level) bispec- 
trum given by equation (|25p breaks down at a rel- 
atively small fc due to non-linearities of gravitational 
growth as well as due to non-linear or even non-local 
bias. Therefore, a model of the bispectrum that takes 
into account higher-order perturbations would be neces- 
sary to push fcmax further. New promising techniques 
such as a renorma lizcd perturbation theory approach 
(ICrocce fc Scoccimarro 2006b. a. 2007t iMcDonaldl l2007t 
iMatarrese fc Pietronil 120071 ) may be used to obtain bet- 
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ter predictions for the power spectrum and bispectrum. 
Further progress is requ ired particularly fo r understand- 
ing redshift distortions ()Scoccimarrcill20d^ . 

In this paper we use a very simple prescription for get- 
ting fcjnax- We choose fcmax SO that a{Raiin, z) = 0.5 

and fcinax = 7r/(2_Rniin). The main motivation for this 
choice being that for small perturbations in the mat- 
ter distribution, say, a{R,z) < 1, one may reasonably 
expect that an analytical model for non-linearities is vi- 
able. Note that /cmax derived in this way depends on z, 
as a{R,z) = a{R,0)D{z)/ D{0). 

An alternative, much more conservative estimate of 
fcmax could be given by requiring that the error on bi de- 
rived from the tree-level bispectrum for some fcmax does 
not exceed the higher-order ("1-loop", or 4th-order per- 
turbation) corrections in perturbation theory to the re- 
duced matter bispectrum at the same /cmax- Specifically, 
one may use 



where ly = 6c/(j{M, z) with 5c = 1.686, and 



Abi ^ AQ^,"'°°P(fcma. 

^1 " Q\T{k^^y) 



(33) 



to determine fcmax- Here, A6i is computed including all 
scales down to fcmax and AQ^~^°°p (fcmax) is the l-loo 
corrections (jScoccimarrol 119971 : IScoccimarro et al]|199 
to the tree-level reduced bispectrum, Q*""®®, evaluated for 
equilateral configurations with fci = fc2 = fca = fcmax- 
This approach, however, makes no use of a large amount 
of information on small scales, and is far from being op- 
timal. Also, the 1-loop correction to the bispectrum, as 
it is the case for the power spectrum, tends to overesti- 
mate the non-linear behaviour measured in simulations, 
thereby making this approach even more conservative 
than necessary. In Sec. 14. 1.21 we shall compare these two 
approaches as a function of volume and number density. 
For the expressions of the 1-loop corrections to th e re- 
duced bispectrum see, e.g.. iBernardeau et al.l ()2002f ). 

3.4. Fiducial values for the galaxy bias parameters 

The galaxy bias parameters, bi and 62, depend on a 
number of factors, including galaxy populations, lumi- 
nosities, and redshifts. On the other hand, the bias of 
dark matter halos, which can be calculated from A^-body 
simulations, is understood relatively well. Therefore, the 
galaxy bias can be calculated from the dark matter halo 
bias, if we assume that galaxies form in dark matter ha- 
los. To do this, one nee ds (at least) the f ollowing informa- 
tion: (i) the halo bias (|Mo et al.lll997t ISheth fc TormenI 
Il999| ) , and (ii) how each halo is populated with galaxies, 
that is, the Halo Occupation Distribution (HOD), {N)m- 

We calculate the galaxy bias parameters from the large 
scales expression 



h ~ — [ dMnh{M,z) b'l{M,z){N)M, (34) 

for i — I and 2, where nfi{M,z) is the mass function 
of dark matter halos of mass M at redshift z, b^{M,z) 
is the halo bias function, and the HOD, {N)m, is the 
mean number of galaxies per halo of a given mass, M . 
We shall use the Sheth fc Tormen's formula for nh{M, z) 
(|Sheth fc Tormenlll999l) : 



n,,(M, z) 



diner 



AP d\uM 



(35) 



V TT '■ 



(36) 



with A ^ 0.322, p = 0.3 and q = 0. 707. The halo 
bias parameters, b^ and 62: are given by ()Mo et al.lll997t 
iScocdmarro et al. 2001b) 



b\{M,z\zf) = l + ei+Ei, 



62^A/,z|z/) = -(ei 



-£■1 ) + £2 + E2 , 



(37) 



where z refers to the redshift of observation, while z/ 
refers to the redshift of formation of halos of mass M , 
and 

qv^ - 1 qv'^ qv^ - 3 
ei = — ; , ^i = - ; , (38) 



l + (g!^2)p' 



Ea 

El 



l + 2p 



2ei, 



(39) 



with Sf = dcD{z)/D{zf), where D{z) is the linear growth 
function. In the left panel of Figure [3] we show the halo 
bias functions, &i(M, z) and 62 (-W, z), as a function of M 
and z, in the approximation that the formation redshift 
equals the observation redshift, z — Zf. We shall always 
assume this throughout the paper. 

As for the HOD , we adopt the form proposed by 
[Tinker et all (|2005h : 



fAT\ 1 , f 



(40) 



for M > Mmin and zero otherwise. The parameter 
Mmin represents the minimum mass above which we find 
a (central) galaxy in the halo, while Mi represents the 
mass above which we can find a second (satellite) galaxy. 
Measuring the HOD parameters for subhalo populations 
from several N-body simula t ions at different redshifts 
and densities, IConrov et al.l (|2006l ) found a correlation 
between Mcut and Mi given by 

logio Mcut = 0.76 fogio Ml + 2.3. (41) 

One also finds from Table 2 in' Conrov et all (|2006( ) that 
Ml /Mmin depends on redshift and density only weakly; 
thus, for simplicity we shall keep this ratio fixed at 
logio(Mi/M„i„) = 1.1, and find Mmin from 



dMnh{M,z){N) 



(42) 



M„ 



for a given rig. 

In the right panel of Figure [3] we show the galaxy bias 
parameters, 61 and 62, from equation (|34p as a function 
of redshift for two values of the mean galaxy density, 
rig = 5 X 10~^ and 5 x lO""* Mpc"^. As expected, 
for a fixed galaxy number density the value of the linear 
bias, bi, increases with redshift. We find that 61 and 62 
are strongly correlated. We shall come back to this point 
in Sec. O 

We admit that these values are derived from very sim- 
plified models without much justification. We need ob- 
servational data to determine the true bias parameters 
for high-z surveys eventually, although we do not have 
sufficient data for doing so yet. Nevertheless, we find 
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Fig. 3. — {Left panel) The halo bias functions, b'l{M, z) (soUd linos) and ^) (dashed lines), as a function of the mass, M, for z = 

(thick lines) and 2 = 1 (thin lines) in the approximation that the formation redshift equals the observation redshift, z = Zf. (Right panel) 
The galaxy bias parameters, bi and b2, for the mean galaxy density of = 5 X 10~^ (continuous lines) and rig = 5 X 10~* h'^ Mpc~^ 
(dashed lines) 



our approach useful for our purpose of deriving the fidu- 
cial values of bi and 62 with a realistic redshift evolu- 
tion, particularly for redshift surveys spanning a wide 
range in redshift, for which one has to consider a set of 
redshift bins and assume different fiducial values for bi 
and 62 at different z. We note that 61 obtained from 
our method agrees with those obtained in the previous 
work by assuming crg.o — 1 and rip :^ 5 x 10 ~"^ Mpc 



(|Seo fc EisensteinilMoa iTakada et al.ll2006D . 

In Sec. [5] we shall show how one can extend this sim- 
ple picture by introducing a redshift dependence in the 
HOD, and how one can make use of the information on 
galaxy bias derived from the bispectrum to constrain the 
HOD parameters directly. 

4. RESULTS 

In this section we present the results from our Fisher 
matrix analysis of the galaxy bispectrum. We first 
study how the derived constraints on the galaxy bias pa- 
rameters and primordial non-Gaussianity depend on the 
choice of fcmax, taking into account the two approaches 
discussed above. We then study how the constraints de- 
pend on the survey volume and redshift. Finally we shall 
apply our method to make forecasts for several current 
and proposed redshift surveys. 

4.1. Dependence on fcmax,. volume, redshift and number 

density 

4.1.1. fcmax 

As mentioned in the previous section, constraints on 
galaxy bias and primordial non-Gaussianity would de- 
pend strongly on fcmax, the smallest scale included in the 
analysis. 

As an example, we consider sample surveys at two red- 
shifts: the median redshifts of (i) z = 1 and (ii) z = 3. 



Each has the volume oiV — 10 Gpc and the number 
density of = 5 x 10~^ /i^ Mpc~^. (The total number 
of galaxies in the survey volume at each redshift is 50 
million galaxies.) The bias parameters are (i) 61 = 1.5 
and 62 = 0.035 at z = 1, and (ii) hi = 2.6 and &2 = 2.1 
at z = 3. 

In the upper panels of Figure |4] we plot the marginal- 
ized, l-cr, fractional errors on hi and 62 at z = 1 and 
3, assuming Gaussian initial conditions, i.e., /nl — 0. 
We observe an interesting effect: a fractional error on bi 
improves at lower z, while that on 62 improves at higher 
z. (Note that this statement is true only when the same 
fcmax is used at both redshifts. See discussion below.) This 
can be understood as follows. Let us recall the form of 
the galaxy reduced bispectrum (Eq. [19] ) : 



,(fcl,fc2,fc3) 



1 62 
— Q(fcl,fc2,fc3) + 73. 

bi b{ 



Now, Q on the right hand side is independent of z at the 
tree-level when initial fluctuations are Gaussian. There- 
fore, the first term falls as I/&1 at higher z where 61 is 
larger (Fig. [3]). On the other hand, the second term ac- 
tually grows as z: for the current example &2/&1 = 0.016 
at z = 1 and 0.31 at z = 3. Therefore, our sensitivity to 
62 grows with z, while our sensitivity to bi declines with 
z. 

Let us study more quantitatively the sensitivity to bi . 
In the limit of linear bias, 62 = 0, a signal-to- noise of the 
reduced bispectrum of equilateral configurations is given 

by 



Ql{k) 



Vb Bl{k,k,k-z) 



Pl{k;z) 



AQ2(fc) [aP{l3)Y kfs^23 

(xD^{z). (43) 
We expect, therefore, that a signal-to-noise of the bis- 
pectrum from gravitational instability declines with z. 
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Fig. 4. — {Upper panels) Predicted errors on galaxy bias parameters vs the maximum wavenumber, femax- The dashed and solid lines 
show the prediction for a galaxy survey at 2: = 1 and 3, respectively. Each survey is assumed to have the survey volume of V = 10 h~'^ Gpc^ 
and the number density of Tig = 5 X 10~^ Mpc""^. The left panel shows the marginalized 1-a errors on the linear bias, bi, while the 
right panel shows the non-linear bias, 62- Both assume Gaussian initial conditions, /jvz, = 0. The vertical lines show fcmax as determined 
from cr{R; z) = 0.5 for each redshift (see Sec. 3.2). (Lower panels) Predicted errors on primordial non-Gaussian parameters vs fcmax. The 
left panel shows the marginalized 1-a errors on the local model, , while the right panel shows the equilateral model, . The bias 
parameters have been marginalized. 



resulting in an increasing error on 61 at higher z. 

In practice, however, we predict that galaxy surveys at 
higher z should result in better determinations of both bi 
and 62- The reason is quite simple: fcmax at higher z must 
be larger than that at lower z. In the upper left panel of 
Figure[4]we show fc,nax as determined from a{R; z) = 0.5: 
fcmax = 0.17 /iMpc"^ at z 1 and fc^ax = 0.47 /iMpc"^ 
at z = 3. The difference is clear: when the modes up to 
fcmax are included, a survey at z = 3 yields an error on 
hi that is a factor of 5 better than that at z = 1. As for 
62, a survey at z = 3 does better by nearly two orders of 
magnitude. 

How about primordial non-Gaussianity? In the lower 
panels of Figure d] we show the predicted errors on 
and /^^, marginalized over 61 and 62- We find that 



the difference between z = 1 and 3 is negligible at the 
same fcmax- This is a consequence of the fact that a 
signal-to-noise for the primordial bispectrum component 
is not, in the first approximation, redshift dependent. 
For equilateral configurations one finds 

Ql{k) [aE{P)Y Vb B]{k,k,k;z) 



oc constant 



Plik;z) 



(44) 

where we considered only the primordial term, Qi. Nev- 
ertheless, we still predict that galaxy surveys at higher z 
should result in better determinations of both and 
/]y\, as fcmax rnust be larger at higher z and therefore 
many more modes are available for the analysis at higher 
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Fig. 5. — Predicted 1-cr errors on galaxy bias and primordial non-Gaussianity vs the survey volume, V, and redshift, z. The short-dashed, 
long-dashed and solid lines show V = 1, 10 and 100 Gpc'^, respectively, with the galaxy number density of = 5 X 10"'^ Mpc~'^. 
{Upper panels) Fractional errors on the linear bias, bi (left), and non-linear bias, 62 (right), for Gaussian initial conditions, /jvi = 0. 
{Lower panels) Errors on primordial non-Gaussian parameters, (left) and f'^j^ (right), marginalized over bi and 62- 



Results from a much more conservative estimate of 
fcmax (eq. [33| ) will be given near the end of the next 
section. 

4.1.2. Volume, redshift, and number density of galaxies 

Dependence of the predicted errors on volume is 
straightforward: it depends simply on 1/W. Depen- 
dence on 2; is a combination of two effects: (i) how a 
signal-to-noise for a given fcmax grows with z, and (ii) 
how fcmax grows with z. Finally, the number density 
of galaxies determines a signal-to-noise on small scales, 
where the shot noise plays an important role. In partic- 
ular, very high-z surveys at, e.g., z > 3, do not add very 
much if the number density of galaxies is too low. 

In Figure [5] we show how the predicted constraints on 
the galaxy bias parameters, bi and 62, improve with the 
survey volume, as a function of the median redshift, z. 



We used Ug = 5x 10~^ Mpc""^ for the number density 
of galaxies. The fiducial values of bi and &2 are calcu- 
lated for each z from Figure [3] We used fcmax determined 
from a{R, z) = 0.5 for a given z. Since fcmax grows as z 
increases, the predicted constraints on bi and 62 also im- 
prove as z increases. A spike dX z ^ 0.8 in A62/&2 is 
a numerical artifact of &2 being very close to zero. The 
dependence on volume is given simply by \/^/V. 

In the lower panels of Figure [5] we show the predicted 
constraints on primordial non-Gaussianity, marginalized 
over 61 and 62- We find that a survey of the size 
y - 1/1-3 Qp(,3 at z ~ 4 - 6 or y ~ 10 Gpc^ at 
z ^ 1 — 2 is as sensitive to as the CMB data from 
Planck. A more ambitious design, e.g., V ^ IQ Gpc'^ 
at z > 2, can achieve A/]^£- ~ 1, although it depends on 
the number density quite strongly. An even more am- 
bitious design, V ~ 100 ft,""^ Gpc'^, would enable us to 
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Same as Figure [5] but for a smaller number density of galaxies, = 5 X 10 * Mpc 



detect the primordial bispectrum from ubiquitous non- 
Gaussianity "floor" fro m the second-order evolution of 
primordial fluctuations (jBartolo et al.ll20()5h . 

Constraints on the equilateral type of non-Gaussianity 
suffer from a stronger degeneracy between the primordial 
bispectrum and the non-linear gravitational evolution as 
well as non-linear bias, and thus the predicted errors on 
are an or der of magnitude larger than those on . 
(See Sec. 14.21 for more detail.) Nevertheless, a survey of 
V ~ 10 h^'^ Gpc'^ at z ~ 1 should provide a constraint 
that is comparable to that from the WMAP 3-yr data. 

We also computed the Fisher matrix for an all-sky 
survey from z — to 5. We divided the entire red- 
shift range in bins of the size Az = 0.5, and used 
rig = 5 X 10~^ /i^ Mpc~^. We flnd that such a survey 
should provide Afj^l ~ 0.2 and Af^^ - 2. These val- 
ues probably represent the best limits on /nl one can 
ever hope to achieve from galaxy surveys. 

How about the number density of galaxies? When the 



number density is low, the shot noise completely dom- 
inates at small scales, and thus one fails to improve a 
signal-to-noise by increasing femax- This suggests that 
very high-z galaxy surveys do not add much if the num- 
ber density of galaxies is too low. In Figure [5] we show 
the case for Ug = 5 x lO-^/i^Mpc"^. Clearly, our sensi- 
tivity to all of bi, &2, /]vL ^-iid does not improve at 
all beyond z ~ 3. Therefore, it makes sense to conduct 
very high-z surveys, only if one can detect more than 
Hg ~ 10-3/i3Mpc"^. 

How robust are these results? The most uncertain pa- 
rameter in our analysis is fcmax- What if fcmax is signif- 
icantly lower than that from a{R, z) — 0.5? To address 
this question, we have repeated our analysis using a much 
more conservative estimate of fcmax given by equation 
(j33p . This estimate was derived by throwing away any 
information beyond fcmax at which the tree-level bispec- 
trum becomes inaccurate. This is a conservative estimate 
because we can certainly improve our theoretical predic- 
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Fig. 7. — Same as fig. [5]but for a much more conservative estimate of fcmax given by equation I I33I I. 



tion by going to the higher order, "1-loop" (4th order ) 
calculations (|Scoccimarroll997l : IScocciniarro et al.lll998D . 
We show the results in Figure [71 We find significantly 
weaker constraints; for example, fractional errors on bi 
from a survey of F = 10 h^^ Gpc^ now go from 9% to 
a few percent from z = 0toz = 5, ora factor of ~ 8 
and 20 weaker constraints at z ~ 1 and 5, respectively. 
For this choice of fcmax the shot-noise is not very impor- 
tant because we are restricted to a fairly small k already. 
Therefore we obtain similar results for a lower density, 
rig = 5 X 10-'^ Mpc~^ . 

4.2. Parameter degeneracy 

Are galaxy bias and primordial non-Gaussianity in- 
dependent? In Figure [5] we show the 2-d joint con- 
straints (95% C.L.) on (6i, 62), {bi, Inl), and (62, 
Jnl), marginalized over Jnl, ^2, and 61, respectively. 
The survey parameters are = 10 Gpc^ , rig — 



5 X 10-3 Mpc" 



1 (left panels) and 3 (right pan- 



els). The fiducial values of bias parameters are b\ = 1.5 
and &2 = 0.035 at z = 1, and 61 = 2.6 and &2 = 2.1 at 
z = 3. 

We find that is not degenerate with 61 or &2, which 
is a very good news; however, reveals a rather strong 
degeneracy with both &i and 62- Therefore, the equilat- 
eral model turns out to be much harder to constrain by 
CMB or galaxy surveys. 



4.3. Current and proposed redshift surveys 

We are now in a position to apply our Fisher ma- 
trix analysis tools to several current and future redshift 
galaxy surveys, both at low and high redshifts. For the 
sake of simplicity, and to allow for easier comparison, we 
shall assume that each survey is characterized uniquely 
by its survey volume, V , redshift range, and galaxy num- 
ber density, Ug. In other words, we shall ignore com- 
plications related to the specific geometry and selection 
functions. 
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Fig. 8. — Two-dimensional joint 95% C.L. constraints on galaxy bias and primordial non-Gaussianity from 2 = 1 (left) and z = 3 (right). 
The top left, bottom left, and bottom right panel show the joint constraints on (fei, 62)1 (^li /jvl)> and (621 /jvi)i marginalized over /atz,, 
62, and bi, respectively. The assumed survey volume is y = 10 Gpc^, while the number density is rig = 5 X 10-3 h^Mpc-^. 



A galaxy survey of a large volume and a relatively low 
galaxy density (but not too low — there should be at 
least Ug ~ 1/P(fcmax) galaxies for the optimal survey 
efficiency) is necessary to detect the baryon acoustic os- 
cillations (BAO) in the power spectrum, which may be 
used to constrain the nature of dark energy. We point 
out that such a survey should also provide competitive 
constraints on primordial non-Gaussianity. 

We shall consider the following survey designs. Table [1] 
tabulates V, Ug, z, fcmax, bi and &2 of these surveys. We 
calculate fcmax from a{R,z) = 0.5 and fcmax — 7r/(2i?) 
(Sec. [3]), and 61 and &2 from a halo approach given in 
Sec. [H 

(i) Low z surveys include two on-going surveys and 
one planned survey: 

— The Sloan Digital Sky Survey (SDSS) main 
sample (for which a detailed analysis of the 
expected constraints on galaxy bias and pri- 
mordial non-Gaus s ianity has been given in 
IScoccimarro et al.l (|2004f )) at z = and V = 
0.3 h-^Gpc^. 

— The SDSS Luminous Red Galaxy (LRG) sam- 
ple at z = 0.35 and V = 0.72 h-^ Gpc^. 

— A proposed extension of SDSS, APO-LSS sur- 
vey at z = 0.35 and V = 3.8 h-^ Gpc^ 

(ii) Intermediate z surveys include three planned sur- 
veys: 

— The Hobby-Eberl y Dark Energ y Experi- 
ment (HETDEX) (Hill et al.ll200l . targeting 
Lyman-a emitters (LAE) between z = 2 and 
z = A. We assume V = 2.7/i"^Gpc^ and 
Ug = 5x 10-'' /iHl^clf^following the "G2" 
design given in iTakada et al.i (|2006. ) for com- 
parison. 



— Two redshift surveys with the planned Wide- 
Field Multi-Object Spectrograph (WFMOS), 
one detecting 2 million galaxies at 0.5 < 
z < 1.3 in 2,000 deg^ (WFMOSl; V = 
4/1-3 Gpc^), and the other detecting 600,000 
galaxies at 2.3 < z < 3.3 i n 300 deg^ (WF- 
M0S2; V = 1/1-3 Gpc^) (jGlazebrook et al.l 
I2005| ). To facilitate comparison we assume 
the same number density of galaxies for both, 
rig = 5 X lO-^'/i^Mpc-^. 

— The Advanced Dark Energy Physics Telescope 
(ADEPT) mission, a space-based redshift sur- 
vey of 100 million galaxies at 1 < z < 2 in 
28,600 deg2. We assume V = lOO/i^^Gpc^ 
and Ug = lO-^/i^Mpc'^. 

(iii) High z surveys are represented by the Cosmic In- 
flation Probe (CIP) mission, a space-based red- 
shift survey targeting Ha emitters at 3.5 < z < 
6.5 (iMelnick et all I2004D . We assume V = 
3.4 h-^ Gpc^ and Ug = 5 x 1 0'^ Mpc~3 follow- 
ing the "SG" design given in ITakada et all (|2006[ ). 

The 8th and 9th column in Table [1] tabulate A6i and 
Ab2 for Gaussian initial conditions {/nl = 0), the 10th, 
11th, and 12th column tabulate Abi, A62, and A/]^£-, 
and the 13th, 14th, and 15th column tabulate A&i, A62, 
and Af^^. For HETDEX, WFMOS, ADEPT, and CIP 
we provide A61, A62, ^/jvi'i ^I'nl from each red- 
shift bin as well as A/]^£- and A/^^ from a combined 
analysis of all bins. 

The predicted constraints on galaxy bias range from 
a few percent for APO-LSS, WFMOS and HETDEX to 
less than 1% for ADEPT, the latter being significantly 
better owing obviously to a larger survey volume. As 
we have mentioned already in Sec. 14. 2[ inclusion of 
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TABLE 1 

Survey parameters (V in units of h'^Gpc'^, Ug in units of h^Mpc"'', z), maximum wavenumber in the analysis (fcmax in units 

OF /iMpc"-*-), FIDUCIAL VALUES OF GALAXY BIAS (fel AND 62), AND MARGINALIZED 1-(T CONSTRAINTS FROM THE FiSHER MATRIX ANALYSIS 

OF THE REDUCED BISPECTRUM IN REDSHIFT SPACE FOR (Tg = 0.9 AND Tls = 1. 





V 






fcmax 




bo 


A61 


Ab2 


A61 


A 62 


A J-loc. 


Abi 


Af)2 




SDSS 


0.3 


30 





0.09 


1.19 


—0.10 


0.270 


0.151 


0.309 


0.151 


255.5 


0.450 


0.421 


1775 


LRG 


0.72 


1 


0.35 


0.11 


2.14 


0.96 


0.209 


0.348 


0.223 


0.353 


113.4 


0.338 


0.726 


998 


APO-LSS 


3.8 


4 


0.35 


0.11 


1.69 


0.21 


0.069 


0.068 


0.071 


0.069 


34.9 


0.108 


0.160 


386 


WFMOSl 


1.6 


5 


0.7 


0.14 


1.87 


0.45 


0.076 


0.096 


0.080 


0.096 


41.0 


0.123 


0.216 


435 




2.4 


5 


1.1 


0.18 


2.16 


1.00 


0.047 


0.081 


0.048 


0.081 


23.1 


0.076 


0.175 


266 




combined 










" 


" 


" 


" 


20.1 


" 


" 


227 


ADEPT 


45 


1 


1.25 


0.20 


2.97 


3.44 


0.020 


0.063 


0.021 


0.063 


6.1 


0.031 


0.111 


73 




55 


1 


1.75 


0.26 


3.44 


5.43 


0.017 


0.066 


0.017 


0.067 


4.5 


0.025 


0.112 


53 




combined 










" 


" 


" 


" 


3.6 


" 


" 


43 


Wr MUbz 


0.5 


5 


2.55 


0.38 


3.27 


4.64 


0.058 


0.220 


0.060 


0.223 


25.7 


0.094 


0.406 


256 




0.5 


5 


3.05 


0.48 


3.64 


6.39 


0.056 


0.253 


0.058 


0.255 


22.1 


0.087 


0.439 


215 




combined 


















16.8 






164 


HETDEX 


0.68 


5 


2.25 


0.34 


3.05 


3.70 


0.051 


0.172 


0.053 


0.174 


23.6 


0.083 


0.326 


244 




0.69 


5 


2.75 


0.42 


3.42 


5.32 


0.049 


0.199 


0.050 


0.201 


20.0 


0.077 


0.357 


202 




0.67 


5 


3.25 


0.53 


3.79 


7.16 


0.050 


0.237 


0.051 


0.238 


18.0 


0.076 


0.401 


177 




0.64 


5 


3.75 


0.65 


4.14 


9.20 


0.053 


0.291 


0.054 


0.292 


17.1 


0.079 


0.469 


163 




combined 


















9.6 






95 


GIF 


1.26 


50 


4 


0.71 


3.16 


4.12 


0.010 


0.036 


0.010 


0.037 


4.7 


0.016 


0.066 


51 




1.13 


50 


5 


1.03 


3.72 


6.76 


0.010 


0.047 


0.010 


0.048 


4.0 


0.015 


0.079 


40 




1.02 


50 


6 


1.46 


4.26 


9.90 


0.011 


0.066 


0.012 


0.066 


3.8 


0.016 


0.102 


36 




combined 


















2.4 






24 



does not degrade sensitivity to galaxy bias, whereas 
does degrade it significantly. 

Our prediction for SPS S should not be cornpared di- 
rectly to those derived bv IScoccimarro et al.l ()2004f ). as 
we have ignored covariance between different bispectrum 
configurations due to non-linear effects and survey ge- 
ometry (Sec. I3.2p . which was included in their work. In 
this sense our analysis is more optimistic; however, in 
the other sense our analysis is in fact more realistic than 
theirs. We have used fc^ax ^ 0.09 /iMpc"^ for SDSS, 
which is significantly more conservative than their value, 
fcmax — 0.3/iMpc~^. This has made our predicted errors 
much weaker than theirs: they obtained Abi/bi = 0.04 
(ours 0.26) and Af°j; = 145 (ours 256). 

How would these results depend on erg and rig? A lower 
value for the rms density fluctuations, as = 0.75, as re- 
centl y suggested by the WMAP 3yr data (jSpergel et al.l 
[2006h . increases fcmax because non-linearity is weaker for 
a lower erg, and thus improves the constraints. On the 
other hand, it follows from equation that the signal- 
to-noise for a given triangle configuration is proportional 
to (t|, and thus a lower cg results in worse constraints. 
When combined, these two effects result in 20% and 
30% stronger constraints on galaxy bias and primordial 
non-Gaussianity, respectively, for low-redshift surveys, 
whereas these effects cancel for intermediate and high 
z surveys. A departure from a scale invariant spectrum 
has a smaller effect. For Ug = 0.95 we find only 10% 
improvement in bias and non-Gaussianity. 

5. CONSTRAINING THE HOD 

So far we have assumed that the linear and quadratic 
bias, bi and 62, are completely free. For surveys spanning 
a large redshift range, therefore, we had to introduce 
multiple redshift bins, and, as a result, we had to have 
an excessive number of free parameters, two times the 
number of redshift bins. 



In this section we attempt to reduce the number of 
free parameters by using a halo approach. The fiducial 
values of bi and &2 were derived from a given form of 
the HOD (Sec. 13. 4p . If we can parametrize the HOD by 
fewer parameters than two times the number of redshift 
bins, the model has more constraining power. 

In the limit that the evolution of bias is given by the 
mass function, we may approximate that the HOD is 
independent of z. We make a minimal extension of the 
HOD that we used in Sec. 13. 4| namely, instead of fixing 
a ratio of Mi and Mmin, we let it free: 

logio Mi^a + logio Mjnzn- (45) 

The fiducial value is a = 1.1, as before. We still assume 
that a relation between Mcut and Mi is given by equation 
([4T|) . We then replace bi and &2 at different redshift bins 
with a single parameter a in the Fisher matrix analysis. 

The expected l-cr errors on a and primordial non- 
Gaussian parameters are given in Table [2] in the third 
to fifth columns. We find that the error on is un- 
affected: this is an expected result because is not 
degenerate with galaxy bias. On the other hand, we find 
a significant, about a factor of two, improvement in f'^i^. 
This is due to the fact that the analysis in terms of the 
HOD is equivalent to introducing a theoretical prior on 
a relation between bi and &2, lifting the degeneracy. In 
other words, while in the previous section we allowed bi 
and 62 to vary independently, we are now making use of 
the fact that the Halo Model predicts them to be strongly 
correlated. 

When we have a survey that covers a wide range in z, 
we may be able to constrain more-than-one parameters 
in the HOD. To see how it works we extend the minimal 
model by introducing one more parameter: 

logoff Ml = a + bz + log^QMmin, (46) 

and we assume 6 = as the fiducial value. An analysis 
based on a single redshift bin would naturally lead to a 
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TABLE 2 

Marginalized 1-a errors on the HOD parameters and primordial non-Gaussianity from the Fisher matrix analysis of the 

REDUCED BISPECTRUM IN REDSHIFT SPACE FOR erg = 0.9 AND rig = 1. THE SURVEY PAR AMETERS AND fcmax ARE THE SAME AS IN TABLE [T] 



Survey 1-parameter HOD 2-parameter HOD 





2 


Aa 


Aa 




Aa 




Aa 


Ab 


Aa 


Ab 


A/'„°^; 


Aa 


Ab 




SDSS 





0.20 


0.26 


198 


0.20 


500 








LRG 





0.10 


0.15 


112 


0.13 


363 








APO-LSS 


0.25 


0.064 


0.077 


35.1 


0.066 


120 








WFMOSl 


0.7 


0.050 


0.066 


39.0 


0.057 


134 


540 


780 


570 


820 


41 


1800 


2500 


435 




1.1 


0.024 


0.031 


22.7 


0.029 


91 


1050 


960 


1070 


970 


23 


3100 


2800 


266 




comb. 


0.024 


0.028 


19.6 


0.026 


75 


0.14 


0.14 


0.15 


0.14 


20 


0.15 


0.14 


75 


ADEPT 


0.7 


0.0050 


0.0062 


6.1 


0.0078 


35 


56 


45 


56 


45 


6.1 


116 


93 


73 




1.1 


0.0033 


0.0040 


4.5 


0.0050 


28 


1200 


710 


1200 


710 


4.5 


2400 


1300 


53 




comb. 


0.0033 


0.0034 


3.6 


0.0042 


22 


0.019 


0.012 


0.020 


0.012 


3.6 


0.020 


0.012 


22 


WFMOS2 


0.7 


0.0147 


0.020 


26 


0.021 


120 


83 


33 


83 


33 


26 


177 


70 


256 




1.1 


0.0120 


0.016 


22 


0.017 


110 


789 


259 


789 


259 


22 


1539 


505 


215 




comb. 


0.0120 


0.013 


17 


0.013 


81 


0.11 


0.038 


0.11 


0.038 


17 


0.11 


0.038 


81 


HETDEX 


2.25 


0.015 


0.020 


24.1 


0.021 


111 


3500 


1600 


3500 


1600 


24.1 


7900 


3500 


250 




2.75 


0.011 


0.015 


19.9 


0.016 


99 


1400 


500 


1400 


500 


19.9 


2800 


1000 


201 




3.25 


0.010 


0.013 


18.0 


0.014 


93 


3200 


990 


3200 


990 


18.0 


6100 


1900 


177 




3.75 


0.010 


0.013 


17.2 


0.013 


92 


730 


200 


730 


200 


17.2 


1300 


350 


164 




comb. 


0.005 


0.007 


9.7 


0.008 


49 


0.034 


0.11 


0.034 


0.011 


9.7 


0.035 


0.011 


49 


CIP 


4 


0.0032 


0.0041 


4.7 


0.0042 


23 


11.1 


2.8 


11.1 


2.8 


4.7 


24 


6.0 


51 




5 


0.0026 


0.0034 


4.0 


0.0034 


21 


6.9 


1.4 


6.9 


1.4 


4.0 


13 


2.6 


40 




6 


0.0025 


0.0033 


3.7 


0.0034 


21 


13.1 


2.2 


13.1 


2.2 


3.8 


22 


3.7 


36 




comb. 


0.0016 


0.0021 


2.3 


0.0021 


13 


0.010 


0.0020 


0.011 


0.0020 


2.4 


0.011 


0.0020 


13 



large degeneracy between a and b. However, one can lift 
this degeneracy by including multiple redshift bins. 

In Table [21 columns 6th to 13th, we present the ex- 
pected 1-a errors on a and b as well as on primordial 
non-Gaussianity. One can clearly see that a and b are 
degenerate within a single redshift bin: combined errors 
are orders of magnitude smaller than those from single 
bins. The determination of primordial non-Gaussianity 
is largely unaffected by the extra HOD parameter. 

6. CONCLUSIONS 

The quest to understand the nature of dark energy has 
recently provided a further motivation for future large 
redshift surveys. It is certain that the study of higher 
order correlation functions of galaxies will be required 
in order to extract maximum cosmological information 
from such large data sets. For instance, as it has been 
recognized for more than a decade, the bispectrum can 
be used to measure non-linearities in the galaxy-mass 
relation. Non-linearity in galaxy bias must be under- 
stood better to meet the high accuracy required for pre- 
cision measurements of the baryon acoustic oscillations 
and their interpretation as a standard ruler, particularly 
for highly biased tracers such as the luminous red galax- 
ies and Lyman-a emitters. 

We have shown that, with the most conservative as- 
sumption about the maximum wavenumber, fc,„ax, used 
in the analysis, the bispectrum measured from a galaxy 
redshift survey should yield a fractional error on the 

linear bias of order 0-l\J h^^ Gpc^ /V at z = to 

0.05^ h-^Gpc^/V at z ^ 6. This is an extremely con- 
servative limit, however, as it assumes no understanding 
of even the mildest non-lincaritics in the dark matter and 
galaxy biasing evolution. 

An analysis that includes all configurations down to 
mildly non-linear scales, within which a{R) < 0.5 is sat- 



isfied, should yield more than an order of magnitude bet- 
ter determination of bias parameters, both linear and 
non-linear, owing to the large number of configurations 
available at smaller scales. This is precisely where in- 
termediate to high-z galaxy surveys play a leading role: 
non-linearity is much weaker at higher z, and therefore 
we can access to a large number of modes on small scales. 

While we can study non-linearities in the matter distri- 
bution by means of new techniques based on perturbation 
theory and of N-body simulations, a simple, local descrip- 
tion of galaxy bias may have to be improved further. In 
this perspective, our results show that a large amount 
of information can be extracted from higher-order cor- 
relations which, in turn, may be used to constrain more 
sophisticated models of galaxy biasing. 

We have also shown that the bispectrum from large- 
volume, high-redshift surveys is highly sensitive to pri- 
mordial non-Gaussianity. The CMB observations have 
been the best probe of the Gaussian nature of primor- 
dial perturbations so far, and the Planck satellite would 
be quite close to the ideal experiment limit. On the other 
hand, a redshift survey of the large scale structure actu- 
ally contains much more information than CMB, as the 
number of modes available from the three-dimensional 
fluctuations is vastly larger than that from the two- 
dimensional temperature and polarization anisotropics. 

Not only can they provide independent constraints on 
scales smaller than those probed by CMB, but also their 
constraints can be comparable to, if not better than, 
those from CMB. The best limit one can achieve from an 
all-sky survey up to redshift ^ 5 should reach ^0.2 
and '--^ 2, an order of magnitude better than the best 
limits achievable by CMB. The planned surveys such as 
HETDEX and ADEPT should reach the constraints that 
are comparable to those from the WMAP and Planck 
CMB experiments, respectively. It should also be un- 
derstood that galaxy surveys provide the best limits on 
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small scales that are not accessible by CMB. This is par- 
ticularly important when probing scale-dependent non- 
Gaussian models. 

Finally, we have shown that galaxy bias parameters 
modeled by the halo occupation distribution should be 
very useful for lifting parameter degeneracies between 
non-linear gravity, galaxy bias, and primordial non- 
Gaussianity. 

We believe that, while a lot of work still has to be done 
in order to achieve a satisfactory description of the evo- 
lution of non-linearities and galaxy bias at small scales, 

higher order correlation functions will play a crucial role 
in the analysis of future redshift surveys, providing in- 



dispensable information on galaxy bias as well as on the 
nature of primordial perturbations. 
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